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§ 1. Introduction 

Thermodynamic quantities generally depend on many details of the system, and are there- 
fore functions of a large number of variables, however, in the critical regime this dependence 
drops to a smaller number. The resulting functions are referred to as scaling functions. 
Scaling functions generically describe a crossover, wherein the effective degrees of freedom 
of a system can change dramatically as a function of scale. Calculating such scaling func- 
tions in critical phenomena is generally accepted to be much more difficult than calculating 
critical exponents. From a renormalization group (RG) point of view one can think of this 
as being due to the fact that to calculate a critical exponent one only needs a local RG 
linearized around the fixed point of interest whereas, generally speaking, to calculate a 
scaling function one needs a global, non-linear RG that is capable of encompassing more 
than one fixed point. One of the chief difficulties in the latter is developing a "uniform" 
approximation scheme that can describe the crossover between two fixed points perturba- 
tively. Conventional small parameters such as e and 1/N might be adequate for certain 
crossovers but not others. 

Crossovers are induced by some asymmetry parameter which often can be fruitfully 
thought of as an "environmental" variable, such as temperature, system size, magnetic 
field etc. The formalism of "environmentally friendly" renormalization [1,2] offers a quite 
general approach to the solution of crossover problems and the calculation of scaling func- 
tions. Given that the key idea behind the notion of a crossover is the qualitatively changing 
nature of the effective degrees of freedom as a function of "scale" and "environment" it 
implements a renormalization which is capable of tracking the evolving effective degrees of 
freedom in a perturbatively controllable manner. However, it is based on reparametriza- 
tion invariance, as in the original field theoretic RG, rather than Wilson/Kadanoff coarse 
graining. 
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The basic idea is that the relation of the bare couplings to the renormalized ones, which 
can be used to describe, parametrically, a physical system, can be thought of as a coordinate 
transformation in the space of these couplings. In thinking of the renormalized couplings 
as new "coordinates" the conventional field theoretic RG simply expresses the invariance 
of physical quantities under changes of coordinate system. This coordinate invariance is 
an exact invariance of field theory. When calculating a physical quantity perturbatively, 
in spite of the fact that physics doesn't depend on coordinates, the particular choice of 
coordinates can be quite crucial in obtaining a reliable approximation scheme. One can 
understand this clearly in the context of the crossover studied in this paper — dimensional 
crossover induced by finite size effects. 

Reverting for the moment to a coarse graining RG, if we thought of possible coarse 
grainings in a d dimensional ferromagnetic film of size L, one would find that block spins 
of size ^ « I were d dimensional, whilst those of size £ 3> L were d — 1 dimensional. 
Thus this coarse graining procedure reflects a crucial property of the "environment" of the 
system — that it is finite in one dimension. Block spinning is therefore an environmentally 
friendly form of renormalization. In the context of reparametrization, an environmentally 
friendly renormalization is one that yields a set of parameters that give a perturbatively 
reliable description of the crossover. In the finite size context a necessary condition for the 
reparametrization to be environmentally friendly is that it be L dependent. 

In previous papers environmentally friendly renormalization has been used to describe 
various physical quantities for particular crossovers of interest [3 — 5]. In this paper we 
consider dimensional crossover of the specific heat as it is one of the more readily accessible 
quantities from an experimental point of view. In the context of films and experimental 
tests of finite size scaling this was the first experimentally measured quantity [6] . 

From a theoretical point of view the case of a totally finite geometry has been success- 
fully investigated numerically [7] and analytically for both periodic boundary conditions 
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[8] and Dirichlet boundary conditions [9]. This case has recently [10] been further devel- 
oped with particular emphasis paid to the problems created by the existence, for an n 
component order parameter, of massless spin- waves (Goldstone modes). In the case of film 
geometries some progress has been made [11] but no theoretical work has been able to 
access the complete crossover other than for two dimensional films [12]. 

The format of the paper is as follows: in section 2 we analyze the connection between 
the specific heat and the vertex functions of a Landau-Ginzburg- Wilson effective Hamil- 
tonian. By choosing as mass parameter tg = A 2 ^ T rjF L ^ we include all the non-analytic 
dependence of the specific heat in the vertex function T^'^\ We then present a renormal- 
ization of the theory as a reparametrization through normalization conditions on certain 
renormalized vertex functions using a fiducial correlation length as our RG scale. Section 
3 is devoted to perturbative calculations. In particular we calculate the specific heat and 
a specific heat effective exponent to one loop. In section 4, by matching to the known 
asymptotic exponents for a three dimensional Ising film, we access the crossover between 
power law behaviour at the three dimensional end and logarithmic behaviour at the two 
dimensional end. The results in sections 3 and 4 are illustrated in the figures. The paper 
ends with our conclusions. 

§ 2. Renormalization of the Specific Heat 

We consider an O(N) symmetric order parameter described by the "microscopic" Landau- 
Ginzburg- Wilson Hamiltonian 

H[<f B ] = £ J d d x Q(V^) 2 + + \tB(x)<PB + ^<Pb ~ H B (x)ip B ^ (2.1) 

which describes a d dimensional film geometry of thickness L. The variable tg when taken 
to be homogeneous has analytic dependence on temperature T, and we choose its origin to 
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be the critical temperature of the film. Hence m? B is determined by the difference between 
the L dependent critical temperature and the mean field critical temperature, i.e. the 
temperature at which the potential in (2.1) acquires a non-zero minimum. \g is assumed 
to be temperature independent. The subscript B refers to bare parameters as distinct from 
renormalized parameters which will be introduced below. We will restrict attention to the 
case when the film also exhibits a phase transition and consider 3<d<4foriV = l, 
and 3 < d < 4 for N > 1. We will present results for periodic, antiperiodic and Dirichlet 
boundary conditions. Note that in the case of periodic boundary conditions our results for 
N = 1 could equally well be reinterpreted to describe the quantum to classical crossover of 
an Ising model in a transverse magnetic field T, where now = V — T C (T) and L = h/T 
[13]. However, in this paper we restrict our considerations to the film geometry. 

The partition function for the the model (2.1) is given by the path integral 



[dip B ]e- n ^-T^ . (2.2) 

The free energy density is F = — y In Z = F — y In Z^qw, where V is the volume and F° 
is the background free energy density obtained after coarse graining from the underlying 
microscopic degrees of freedom to those of the effective field theory description in terms 
of the LGW Hamiltonian (2.1), Z^q^ being the partition function of this Hamiltonian. 
F b is assumed to be an analytic function of the thermodynamic variables. The internal 
energy density is 

dF 

U=F-T— 
dT 

and the specific heat, by definition dU /dT, is given by 

r) 2 F 

C =- T2 df^- ( 2 - 3 ) 

The assumption in working with this LGW Hamiltonian is that the only one of its param- 
eters to retain a dependence on temperature is the mass parameter t B . Thus the internal 
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energy density 

V = rf-^j*x^UP» (2.4) 

and the specific heat 

dt B (y) 
dT 
(2.5) 



C = C»-^j d^T^)G^ x)+ T w J d d x J dS^G^y)- 



where 

G^\x) = (^(x)) and G^\x,y) = ^\x)^\y)) - (<p 2 (x)) (^(y)). (2.6) 

Concentrating on T > T c , where < (p >= 0, and denoting 

Y[t B ] = -\nZ LGW (2.7) 

we have that 

r(°' 1 )(x) = \g^ X \x) and T i0 > 2 \x,y) = -^-G^ 2 \x,y). (2.8) 
So for homogeneous t B we find 

u = u i_ T^rf 1 " 1 (2.9) 

where r^' 1 ^ and T^'^ are to be evaluated at zero external momentum. 

If we wish to incorporate all of the non-analytic dependence of the internal energy and 
the specific heat into iS ' 1 ) and ]?( ' 2 ) respectively, then a natural choice of the dependence 
of t B is 

t B = A 2 (2.11) 

where A is a microscopic mass scale. In the vicinity of the critical temperature the results 
with this variable will be the same as those obtained with the linear measure A 2 ^M. 
With the choice (2.11) the internal energy density becomes 

U = U b - A 2 T c r2' 1} (2.12) 



and the specific heat is given by 



C = C b -^T™. (2.13) 



For an O(N) model G^'^ is manifestly positive and either diverges or goes to zero at the 
critical temperature according to the value of N. Thus we anticipate that r^' 2 -* should 

diverge to — oo or vanish at the critical temperature. Our problem is therefore to calculate 

r (0,2) 
1 B ■ 

The correlation length in the transverse dimensions, = m _1 , we define via the 
second moment of the two point function, On Fourier transforming is obtained 

from 

(2) 

2 T y B '(p,t B (m),X B ,L) 
m = — (2.14) 



d p 2T^\p, //>>(/// ). Xb, L) p2=() 
where p is the transverse momentum and t B {m) is that bare mass parameter which pro- 
duces the inverse correlation length m. The origin for the variable t B (m) is specified by 
requiring that 

rg ) (0,0,A B ,L)=0 (2.15) 

which insures that tg is proportional to T— T C (L) as the critical temperature is approached. 
Changing the mass parameter tg, by changing the temperature in (2.11), allows us to tune 
the correlation length. Note that the physical correlation length of the film geometry (2.14) 
depends on L and will be infinite at the film critical temperature T C (L). 

We will define renormalized parameters by 

t(m, k) = Z^ 2 1 (K)t^(m) and X(k) = Z\{k)\q (2.16) 
and renormalized vertex functions by 

r^ M )(m, k) = zf( K )Z^(K)T^ M \m) + 5 m 5 Mn A^{n) n = 0, 1, 2 (2.17) 



which is just a reparametrization of the original theory, where k is an arbitrary renormal- 
ization scale. 

Contrary to the renormalization of other vertex functions, like and T^, the vertex 
functions r^'^ (n = 0, 1, 2) have to be renormalized additively via 

r (0,n) = Z ™ 2 rf n) + (2.18) 

T and T^ ' 1 ) determine the Gibbs free energy density and the energy density of the LGW 
Hamiltonian respectively * . 

Before discussing the renormalization of r(°' n ) for n = 0,1,2 we will specify the Z J s 
associated with the reparametrization (2.16) and (2.17). Here we will restrict ourselves 
to T > T C (L). The case of T < T C (L) will be considered in conjunction with crossover 
amplitude ratios elsewhere. For T > T C (L) the conditions which specify our Z : s are 



z- 1 = dj?r%(p,t B ( K ),\ B ,L) 



,(2) 



p 2 =0 



(2.19) 



-i _ (p,*b(«),a b ,l) 



(2) 

d p 2T B J (p,t B (K),X B ,L) 
r ( R 4) (0,t B ( K ),As,L) 



p 2 =0 



where the relation between t B (n) and k is specified by 



K 



T^(p,t B ( K ),X B ,L) 



(2) 

d P 2T B (P,t B {K), X B ,L) 

and the origin of t B is fixed by (2.15). Note that the Z : s are obtained from the vertex 
functions of the system specified at an arbitrary, fiducial, transverse correlation length k -1 , 
as opposed to the correlation length of interest, m _1 . Furthermore the conditions are all 



p 2 =0 



(2.20) 



(2.21) 



(2.22) 



For homogeneous t B , H B = and T > T c we use the convention r = y In Z^qw 
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L dependent. As has been emphasized on previous occasions [1,2] such "environmentally 
friendly" conditions are essential in order to obtain a perturbatively controllable description 
of the finite size crossover. 

We define the Wilson functions as the logarithmic derivatives 

1 dZm 

> = V*"* < 2 ' 23 > 

1 dZ, n 2 

v - -z*t (2 - 24) 

The Wilson functions 7^2, 7^ and 7^ are explicitly L dependent and interpolate between 
those of a d and d — 1 dimensional O(N) model in the limits kL — > 00, k — > and 
kL — > 0, k — ► respectively. 

The invariance of the bare vertex functions, r^'^, under the one parameter group 
of reparametrizations indexed by the arbitrary renormalization scale k (they don't know 
which reference correlation length k -1 will be picked to define the reparametrization) yields 
the RG equation 

jy(N,L) AT 

+ (L^2 - -1^ {N ' L) = S m S Ln B^ (2.26) 

where n = 0, 1,2. The equation is inhomogeneous for the three vertex functions T, iS ' 1 ) 
and r( ' 2 ), where the "source" term 

1 ^ dA^ t \ 

B (n) = + nl A (n) (2 _ 27) 

da v 

is finite order by order in the loop expansion. 

The relationship between temperature and n can be obtained by using 

ft 

r( 2 )(t) = / Y^ l \t')dt' (2.28) 

Jo 



and conditions (2.19-2.22) with the definitions of the Wilson functions (2.23-2.25) to find 

t(m, k) = k 2 ^ —(2 - 7 ^) e r( 2 -V)f . (2.29) 
JO x 



We see that 



dt(m, k) , 
« \ K =l^t{m,K). (2.30) 



An important feature of the above is that the determination of 1?( 2 ) by integrating H 2 ' 1 ) 
allows us to bypass the need to determine m 2 B perturbatively. 

In terms of the renormalized vertex functions, the conditions (2.19-2.22) are equivalent 

to 

a^r( 2 )(p, *(«,«), a, =i (2.3i) 

T^(0, t(n, k), A, L,k) = 1 (2.32) 

r^(0, t(n, k), A, L, k) = A. (2.33) 

r( 2 )(0,t(K,^), \,L,k) = k 2 . (2.34) 
We could have replaced (2.32) by the condition 

t(K,K)=K 2 . (2.35) 

This condition together with (2.34) determines a multiplicative renormalization of £g, and 
of (p 2 insertions via a renormalization function Zf. The two renormalization functions Zf 
and Z^2 are different, the latter being determined by (2.32). The quantity 7t = — ^fjj? is 
an analog of 7^2, however, the problem with implementing a condition such as (2.35) in 
perturbation theory is that the resulting Zt involves diagrams with massless propagators, 
some of which are strictly infinite even after the introduction of an ultraviolet cutoff. 
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Defining an effective exponent v ef $ = , one finds 



si'*,,* „ a>T( 2 -v) 



^//-(2-7i) ^^j") • ( 2 - 36 ) 

One can also define what we term a floating exponent, v f = (2 — 7< j3 2)~ 1 . As near a fixed 
point 7^ and 7^2 go to constants we can see from (2.36) that 7t — > 7^2, hence both the true 
effective exponent and the floating exponent interpolate between the same two asymptotic 
values. This may not be evident perturbatively. One can think of the floating exponents 
evaluated in environmentally friendly RG improved perturbation theory as approximations 
to the true effective exponents [2] . Another way of thinking about them is from the point 
of view of a redefined temperature variable in the following way: if one defines t' = tf(t) 
and u f = = (2 — 7^2 ) _1 one finds that 

dint (2-7*0 h x K m 1 ' 

Near a fixed point / — > 1 hence t' — > £. 

The solution of the RG equation (2.26) for the specific heat is 

r(°' 2 )(t(m, k), A(k), L, k) =e 2 $* ^VrM^m, /ok), A(p«), L, «p) 



1 x 



(2.38) 



Reparametrization invariance is now manifest in the fact that the left hand side of (2.38) 
is independent of p, the latter being just an arbitrary rescaling of n. 

We will now discuss some possible normalization conditions for r( ' 2 ), thus specifying 
A^\k). One possible choice is the normalization condition 

r(°> 2 \t(K,K), \(k),L,k) = (2.39) 

which is equivalent to 

aW(k) = -Zl 2 T^ 2 \t B {^\ B ,L). (2.40) 
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The advantage of this condition is that all the "physics", in the sense of the effects of 
all fluctuations, is now purely in the inhomogeneous term. The normalization condition 
(2.39), however, is natural as the ]?( ' 2 ) does indeed vanish in the mean field regime, or at 
least goes to a constant which can be chosen to be zero. Neglecting the inhomogeneous 
term r^' 2 ) being zero is then an invariant statement with respect to RG transformations. 

A methodology which avoids some of the pitfalls of additive renormalization is to relate 
]?( ' 2 ) to the correlation function r( ' 3 ), the advantage of this approach being that the latter 
is multiplicatively renormalizable in d < 6. We have the analog of (2.28) 

r(°> 2 \t)-r( > 2 \ti)= r^ 3 \t')dt'. (2.41) 

Ju 

Using the relation 

r(°< 3 )(t(m, k), A(k), L, k) = e 3 fi 'Vr^^m, P k), X( P k), L, K p) (2.42) 

and the relation between the correlation length and the temperature (2.29) one finds 

r(°< 2 )(t(m, «)) = K d - A r -(2 - ^)Ji {2 ^- 4+d ^ f (°< 3 )(x) (2.43) 

J oo x 

where we have normalized r^' 2 ) to vanish in the mean field limit and 

- (n * r(°' 3 )r( 2 ) 3 

r^(m)= (2.44) 

rl 2 ' 1 ) m d 

It is not difficult to show that in fact (2.43) is exactly the same as the expression (2.38) 
obtained from the additive renormalization prescription with the normalization condition 
(2.39) at t(oo,oo). 



§ 3. Perturbative calculations 



We begin this section by analysing the (3 function for the coupling, as we will perturbatively 
expand all other functions in terms of the solution of this equation . In terms of the floating 
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coupling h [1], chosen to be the leading term in the perturbative series for 7^, one finds, 
for p^ = f3(h, z), to one loop 

/3{h,z) = -e(z)h + h 2 . (3.1) 



The function s, in an obvious diagrammatic notation, is 

e(z) = ^p-2, (3.2) 

depends on d and z = pnL but is independent of N. We take the solution of (3.1) 

e Jz o y ' x 

h (z) = fx - dy — (3-3) 

as our perturbation parameter. 

After solving the equation we specify the arbitrary scale p to be p = and relate 
it to temperature via (2.29) whereupon z becomes In (3.3) the initial coupling is 

then taken to be at a "microscopic" scale k. For d < 4 this microscopic scale can be 
sent to infinity and a universal floating coupling, the separatrix solution h(z) = 4z ^jp 
obtained [2] . Of course, if one is interested in corrections to scaling, as is usually the case 
in comparing with experimental data, then k should be left finite and fitted to the data. 

For periodic boundary conditions one finds 

jr 4?rV ^ , ^Vj^ 

e(z) = 5 - d - (7 - d) n= -°° Z —r^- (3-4) 

n=— 00 v ' 

and the separatrix coupling 

(d-7) 

oo a 9 9 o 

E + 

h(z) = (5 - d)^— ~ (3.5) 

E + 

n=— 00 
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For d = 3 the results are particularly simple 



h(z) = 1 + 



sinhz 



2 2 coth(f) 

e(z) = 1 + — ^ 

sinh z + z 

where, of course, we are now restricted to N = 1. 

We present here the corresponding results for Dirichlet and antiperiodic 
conditions. For Dirichlet boundary conditions 

ttW - 1) / ttV - 1) " 9 
^2 I ^2 

e (z) =5-rf-(7-rf)^ 1 



n=l 



z 2 



and for the separatrix coupling 



OO (rf-7) 

J2a + '-^) 2 



h(z) = (5 -<()"" 1 



OO (d-5) 
+ 2 



n=l 

For antiperiodic boundary conditions one finds 



-(:) 5 d (7 d)" 



d — 9 

7r 2 n(n + l) / : 7T 2 n(n + 



g A _7r 2 n(n + l)^ 7 



72= — OO 



Z 2 



and finally the separatrix coupling 



OO (d-7) 

E (1 + ^^) 2 



//(:) = (5 x 



OO (ti-5) ■ 

£ (i + ^§±^) 2 



71= — OO 
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For d = 3 the results are once again very simple. For the Dirichlet case 

9 9 ( V 2 _ tanhjM 

y y ^ ~+~ sinh2y Z y > 

where y = (z 2 — 7r 2 )2. Even though y(z) has a branch point h(z) is analytic in 2. The 
separatrix coupling is 

2 (1 + 2 I/ _ 2tanhy\ 

= d + • < 3 ' 13 > 

The corresponding results for antiperiodic boundary conditions are 

(314) 

y l (sinhy-y) 



and 

h(v) = n -+■ 

sinhy 



2 



%) = (i + ^)(i--frr:)- (3-15) 



The Wilson function 7^2 is given by 

whilst 7^ = to one loop. Two loop Pade resummed expressions for the Wilson functions 

and the floating coupling, for the case of periodic boundary conditions, can be found in 

[2,5]. Substituting any of the above floating couplings into (3.16) yields 7^2 for the three 

different types of boundary condition. As mentioned corrections to scaling can easily be 

included. For example for d = 3 and periodic boundary conditions 

! zsinh(f) 2 1 sinhzp + ^o coth(f ) sinh^ 

" \ z ) = 77 — r . , ,,„,, — 2, ) + -— . (3.17) 

smhz + z h{zQ) ZQsmh{-f) 1 zq smhz + z 

Turning now to r( ' 2 ), up to two loop order and once again in an obvious diagrammatic 
notation (note that we have made the diagrams dimensionless by pulling out an overall 
scale) r^' 2 ^ is given by 

if = -^W- 4 [0 - A b k"<^>0 2 ] (3.18) 
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the two loop graph with "tadpole" having been absorbed into the one loop propagator 
by the replacement of with pn using (2.22). Implementing the normalization condition 
(2.39) one finds that 

5(2) = -2NK 2 0\ n . p (3.19) 

where the subscript denotes that the diagram is evaluated at the normalization point. Thus 
we see that the one and two loop expressions for in terms of renormalized quantities 
are identical. Explicitly to one loop for periodic boundary conditions one finds 

Bi 2)_ iv r(V)(^)^ ~ 4.V ^ ,„ 2m 

L (2tt) 2 n= _ OG ^ « P 
p(0,2) j s t hus f ounc [ ^ substituting (3.19) and (3.16) into (2.38) to obtain 



2LfC (2tt) — 



(3.21) 



n=— oc 



where the arbitrary scale p, as before, is associated directly with the inverse correlation 
length. Thus we calculate the specific heat and other physical quantities directly in terms 
of the finite size correlation length. Equation (2.29) relating ^ to L and t provides a 
parametric representation of physical quantities in terms of t. 

In the limit p — > only the n = term in the sum is important and one finds 

r(o,2) _ iv(iv + 8) rj^)^ d _ 5+2{5 _ d){m 

2(4 — N)Lk (27r) ^i • 1 • J 

In the same limit one finds p — > (J?)^ -1 , z^-i = (2 — (]\f+f)(5 ~~ c 0) _1 being the d—1 
dimensional correlation length exponent. Hence 

r M _ -^ (JV ,tl r(¥) :!" 4 (4)-^ (3.23) 



2(4 — N)Lk (27r) ^i 
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where a d _ 1 = 7 N „ ^ - w+s - is the d — 1 dimensional specific heat exponent. Similarly, 

in the limit Lup — > oo, p — > the sum can be converted to an integral and one finds that 

r( o,2) _> iv(iv + 8)r(^)^ rf -4 + 2(4- rf) (#±i) (3 24) 

2(4 — AT) ( 27r )d/2 p ^ J 



and p — > (4y) d where ^ is the <i dimensional correlation length exponent. In the bulk 
limit 

r (o, 2) ^ iV(iv + 8)r ( V)^ , 

where = 7 , f , „ x - v N+s ' is the d dimensional specific heat exponent. Thus we see 

2-(^|)(4-d) 

that the specific heat crosses over precisely between the expected d and d — 1 dimensional 
asymptotic forms. 

Note that the amplitude of r( ' 2 ) in the above expressions appears to diverge at iV = 4 
this is an artifact of the one loop approximation. What actually happens is that for d 
between two and four there is some value of iV for which a(N, d) = 0, at this value of N 
and d we expect the specific heat to have a logarithmic dependence on t. For N = 1 this 
occurs at d = 2, however, at one loop the value appears to be independent of d and occurs 
at iV = 4, which is the relevant value for d = 4. 

A plot of the specific heat as a function of correlation length is shown in Fig. 1 for a 
three dimensional Ising film with periodic boundary conditions. The effective specific heat 
exponent defined as 

dlnC (*Of<\ 

= ~ ~dh[t (3 - 26) 

is plotted in Fig. 2 for the same model. Note that in this approximation the asymptotic 
two dimensional value of a eff is 0.5 as opposed to the exact value of zero, obtained from 
the solution of the two dimensional Ising model. This is a weakness of the perturbative 
approach which effects the specific heat exponent in a particularly acute manner. In the 
next section by matching to the known asymptotic exponents of the model we investigate 
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the more realistic behaviour. In the case of a four dimensional O(N) film, in the limit 
^- — > oo, ~^ °°? one finds that 



r (0,2)_ ^ ^ + 



— In —~ 

2 Av 



4-JV 



JV+8 

(3.27) 



16tt 2 \4- iV 

in accordance with known results. Fig.'s 3 and 4 show plots of the specific heat and a eff for 
the four dimensional Ising film. Note the presence in the figures of logarithmic tails at the 
four dimensional end as described by (3.27). Fig. 5 shows a comparison of a ef f for a three 
dimensional Ising film with Dirichlet and antiperiodic boundary conditions. Additionally 
the result for the Gaussian model is plotted in Fig. 6 with periodic boundary conditions. 

§ 4. Crossover to Logarithmic Behaviour in a Three Dimensional Ising Film. 

In this section we will consider the crossover between three and two dimensions for an 
Ising model in a way that is capable of accessing the logarithmic behaviour characteristic 
of the two dimensional specific heat. For the two dimensional Ising model a = 2 — vd = 0. 
The consequent logarithmic behaviour of the specific heat is thus due to a competition 
between v and d. For d = 2 the correlation length exponent v = 1, hence a = 0. Now for a 
three dimensional Ising film with periodic boundary conditions, at one loop the crossover 
is governed by the floating coupling h = 1 + gin z h - . This implies a crossover for u e jj 
between 1/6 and 1/3. By far the biggest error involved in evaluating crossover functions 
is associated with the values of the asymptotic exponents themselves. With this in mind 
one is inclined to try to match the scaling function to the asymptotic exponents. This can 
very simply be done in the case at hand by writing h = A + where now the constants 
A and B will be determined by demanding that as z — > 0, u e jj — > 1 and that as L — > oo, 
z — > one finds v e ff ~ * 0.630. The values 1 and 0.630 are the exact two dimensional and 
three dimensional 6-loop Borel resummed [14] exponents respectively. Thus one finds that 
A = 1.238 and B = 1.762. 
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In Fig. 7 we plot a eff as a function of Inz by substituting our ansatz for h into 
(3.21). Note the logarithmic tail as the two dimensional critical region is approached. 
More interestingly, there is a pronounced bump in the curve which is absent in the one- 
loop approximation. This arises due to a competition between the effects of u eff and the 
effective dimensionality d eff [2] . The bump remains even if one uses a completely different 
interpolating function such as h = A+j^, though its amplitude and width vary somewhat. 
In Fig. 8 we plot analogous results for the case of Dirichlet and antiperiodic boundary 
conditions. Once again the bump is clearly present. In the case of Dirichlet conditions 
however there is also a dip before the bump is reached. Based on previous experience of the 
behaviour of effective exponents with Dirichlet boundary conditions [2] this is not totally 
unexpected. In Fig.'s 9 and 10 we have used instead of the universal floating coupling 
the coupling (3.17). There is now a double crossover; firstly between mean field theory 
and the three dimensional asymptotic exponent and then to the asymptotic behaviour of 
the two dimensional exponent. In Fig. 9 we plot the result for the case where we do not 
match to the exact two dimensional exponent and in Fig. 10 the result with matching. 
The asymptotic three dimensional regime would most probably be much narrower than 
that shown. This can be very easily modeled by adjusting the initial condition for the RG 
flow. In the case at hand, we have, for the sake of clarity, and to emphasize the double 
crossover, left it large. It is clear from the figure how the effective exponent would be 
modified as the well developed three dimensional universal regime is narrowed. 

§ 5. Conclusions 

In this paper, using environmentally friendly renormalization, we have treated the finite 
size crossover of the specific heat of an O(N) model in a <i-dimensional film geometry. For 
N > 1 we considered 3 < d < 4, and foriV = l,3<d<4. We derived expressions for the 
specific heat and an effective critical exponent a eff that were completely regular across 
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the entire crossover, the expansion parameter for the perturbative series being the floating 
coupling h. We considered periodic, Dirichlet and antiperiodic boundary conditions. 

For the crossover from three to two dimensions of an Ising film we saw that one loop 
answers in the asymptotic two dimensional regime were quite poor. As is known, generally 
speaking, perturbation theory becomes more unreliable as one goes to lower dimensions. 
For the specific heat the problem is particularly acute as the one loop effective specific heat 
exponent was seen to be monotonically increasing whereas, as we know from the solution 
of the exact two-dimensional Ising model, the two dimensional specific heat exponent 
is strictly less than the three dimensional one. Hence we could say that the one loop 
approximation is failing to capture a qualitative feature of the crossover in this case. 
To circumvent this problem, and in the knowledge that the dominant source of error in 
calculating scaling functions is the uncertainty in the asymptotic critical exponents, we 
took a more pragmatic line by making an ansatz for the floating coupling constant so as 
to be able to asymptotically match the "known" two and three dimensional correlation 
length exponents. By so doing we were able to access in a very simple way the crossover 
between power law and logarithmic behaviour in the asymptotic regime, finding that the 
resultant crossover curve had a very interesting bump. We also analyzed the crossover to 
mean field theory thereby accessing a double crossover governed by three different fixed 
points. Our global, environmentally friendly RG captured all of these fixed points in one 
uniform approximation scheme. 

The crossover between three and two dimensions for N > 1, and in particular for 
N = 2, are potential problems that could be analyzed using the techniques of this paper. 
The latter being a problem of longstanding interest for experiments with liquid helium 
confined to a film geometry [6,15]. We hope to return to these issues in the future. 
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